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Abstract 

The two-dimensional Jordan- Wigner transformation is used to investigate the 
zero temperature spin-Peierls transition for an anisotropic two-dimensional 
XY model in adiabatic limit. The phase diagram between the dimerized (D) 
state and uniform (U) state is shown in the parameter space of dimensionless 
interchain coupling h (= J±/J) and spin-lattice coupling rj. It is found that 
the spin-lattice coupling rj must exceed some critical value rjc in order to reach 
the D phase for any finite h. The dependence of r]c on h is given by — 1/ln/i 
for h ^ and the transition between U and D phase is of first-order for at 
least h > 10"^. 
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The fact that the S = 1/2 Heisenberg spin chain is unstable towards dimerization when 
coupled to an elastic lattice is known as spin-Peierls (SP) transition This occurs because 
dimerization opens a gap in the spin excitation spectrum and lowers the total magnetic 
energy by a greater amount than the increase in elastic energy due to lattice deformation. 
Such a transition was first suggested by analogy with the conventional Peierls transition in 
linear conducting chains and later observed in several organic compounds such as MEM- 
(TCNQ)2 0]. It has attracted renewed attention since the discovery of inorganic compound 
CuGeOg in 1993 

In the adiabatic limit, it is well understood for exactly one-dimensional (ID) XY or 
Heisenberg model that the dimerized (D) state, relative to the uniform (U) one, is always 
stable for arbitrarily weak spin- lattice coupling in the ground state In presence of some 
extra mechanisms these two states may compete leading to interesting phase diagram. Up 
to now there are several mechanisms which were found to destroy the D phase, for example, 
small Ising anisotropy Another mechanism is quantum lattice fluctuation which was 
intensely studied recently Beyond adiabatic approximation it was found that the spin- 
phonon coupling must be larger than some nonzero critical value for the SP instability to 
occur |^J§]. The interchain exchange coupling is also a possible mechanism This 
mechanism is often not negligible for real SP materials, especially for the material CuGeOs 
where a relatively large interchain coupling, about J_l/J = 0.1 (J represents the intrachain 
exchange) was found experimentally [|10|- Theoretically the importance of the interchain 



coupling was also stressed by many authors in description of CuGeOa [pJ]-[T3[|. Therefore 
it is necessary to study the SP transition for a general quasi-lD system with inclusion of 
interchain coupling. This should be discussed by answering a central problem how the 
interchain coupling modifies the dimerized ground state of the exactly ID case. 

In this work we propose to study the above problem through a two-dimensional (2D) 
XY model with anisotropic intrachain and interchain coupling, and leave the generalization 
to the case of the corresponding Heisenberg model for later study. The reason is two-fold: 
i. the problem is simplified; especially the ID limit can be exactly solved, ii. the XY model 



is believed to contain the essential elements for the spin-Peierls transition . We point out 
that the similar problem based on the Heisenberg model was studied previously by some 
authors J^,]^. However, they all treated the interchain coupling within the so called chain 
mean-field theory and actually still studied a ID system. Usually the effect of interchain 
coupling is underestimated in this way. Here, the two-dimensional (2D) Jordan- Wigner 
(JW) transformation |T^-[l7i3 is adopted so that the intrachain and interchain exchanges 
may be treated on an equal footing from the beginning. 

In adiabatic approximation, we begin with the following anisotropic 2D XY model on a 
square lattice: 

where S^j"^ represents the x/y component of the spin-1/2 operator at site fij = ix + jy 
{x, y are unit vectors along the x- and y-axis, respectively), J{J±) is the intra(inter)chain 
exchange coupling along the x{y)-a.xis, a represents spin-lattice coupling constant and K is 
the intrachain elastic constant. Here the lattice displacement is assumed as Uij = {—ly^^u, 
i.e., corresponding to a wave vector (vr, vr), which was found to be the case in CuGeOs by 



neutron-diffraction measurements |I8|. In the following we take J as unit of energy and 
use several dimensionless parameters: the dimerization parameter 6 = au, the interchain 
coupling h = J±/J and the spin-lattice coupling rj = a^J/K. The total number of lattice 
sites is A^. 

The spin operators may be represented as fermions by use of the 2D JW transformation, 
which has the advantage that all spin commutation relations, as well as the so called spin 
on-site exclusion principle are automatically preserved [^. This method has been proved 
well for application to real materials p5|,|9|. For concreteness, the following formulas are 
adopted |16|,|1] 



= Cij exp{i(pij}, = exp{-i(p,j}clj, = cljCij - 1/2 (2) 

with 



i—1 CO j — 1 

n 



1=0 m=0 m=0 

where c|j(cjj) is creation (annihilation) operator for a spinless electron at rij and riij = 
cljCij is the particle number operator. The original Hamiltonian can then be transformed 
into: 

H = Y.[1 + (-l)^+M](e*-c^c,+i,, + h.c.)/2 + h E(cL-Q,,+i + h.c.)/2 + N6y2r^ , (3) 

where the phase factor e^^^-^ = e^'^'^'+i-'"'^''^^ describes an effective gauge field acting on spinless 
fermions. Note that the Hamiltonian (j^) is exact and it may recover to the corresponding 
one derived from ID JW transformation as /i ^ 0. When going on one has to treat the phase 
factor approximately. Following previous work, we select a configuration (see Figs. 1 in Refs. 
14, [T5|JT7| ) where the phase factor e*'^'-^ has alternative values e*'^ and 1, i.e., it varies with 



(—!)*+■' (note that this coincides with the above dimerization pattern). This configuration 
ensures that each elementary plaquette encloses a net flux of vr |14,T^. Then the Hamiltonian 



can be rewritten as follows in terms of fermion operators e and / corresponding to the two 
sublattices A and B respectively (constants irrelevant to S are ignored) 

H=l i: [-(1 - Wll,e^, + h.c.) + (1 + 5)(e^/,+i,, + h.c.) 

A 

+Kfh-iei, + + h.c.)] + N5^/2^ . (4) 



In the following we want to study the SP transition at T = based on the above Hamiltonian. 
Depending on the value 6* which minimizes the ground state energy, the system may be in 
D phase {6* 7^ 0) or U phase {6* = 0) in the parameter space of interchain coupling h and 
spin-lattice coupling rj. 

The Hamiltonian (H) may be exactly diagonalized. The electronic spectrum is written as 



= ± Y (5 cos kx + h cos kyY + sin^ within the (reduced) Brillouin zone: — vr < k,j. ± ky < 
TT, and the ground state energy is simply given by Eqs = — Z^k kk| +N5'^ /2rj. Before contin- 
uing we add a comment here on the above used JW transformation. In the definition of 0jj 
in Eq. (S) the summation is selected over the lattice sites along the direction perpendicular 



to the chains so that the phase factor appears in the intrachain hopping term as shown in 
Eq. One may also think to do the summation along the chain direction, i.e., re-define 
4'i,j = J2iri=o I^/^o ^i,m + J2\Zo ) and then the phase factor will appear in the interchain 
(rather than intrachain) hopping term. With the similar treatment for this phase factor as 



done above one may obtain an alternative spectrum ±-y/ (5 sin kx + h sin ky)'^ + cos^ k^. How- 
ever, the ground state energy Eqs (i.e., the summation over the Brillouin zone) is not changed 
by the new spectrum, nor all the subsequent results derived from it. Similar discussion was 



given in Ref. ||16l. Then we go on to search for the value S* which may be obtained from 



the condition dEcs/dS = 0, yielding the following equation: 

^ .TT .TT-fc, cos kx + h cos ky) cos k^ 

6 = — dkx dky ' . (5) 

zjc Jo Jk^-TT J (5 COS kx + h COS ky)^ + sin kx 



It is interesting to notice that Eq. (||) will reduce to the corresponding rigorous equation 



for the exact ID case [IS] if we set /i = 0, which means. 



i='-rdkx^ (6) 

^ -^0 J 5^ sin^ kx + cos2 kx 



The Eq. (|^) always has a single nonzero solution for 6 as long as ?7 > ||2D|, which predicts a 
dimerized ground state. For h ^ 0, Eq. (^) may give one or two nonzero solutions 6i {i = 1, 2) 
except for the trivial solution 5 = 0. We need to consider the sign of d'^EQs/d6'^\s=Si, as 
well as compare the energies between at 5 = 6i and at end points 6 = 0, 1 to obtain the 
actual 6*. Without details, we give the results shown in Figs. |l] and |[ The Fig. |l] gives 
the D/U phase diagram in the parameter space {h,ri). As expected, the interchain coupling 
tends to suppress the D phase. As long as the interchain coupling h > 0, the spin-lattice 
coupling rj must exceed some critical value rjc to reach the D phase and the value rjc increases 
monotonically with h. Careful analysis shows that the critical value rj^ has the functional 
form — l/ln/i in the region h 0, which increases starting from zero much faster than any 
power law. Moreover, we have found that the transition from U to D phase with increasing 
1] is of first-order for at least h > 10^^ As shown in Fig. ^ we may see how the 



order parameter 6* jumps from zero to finite value with increase of rj for several h values. 



Although the transition is of second-order at the point h = as can be proven from Eq. 
we speculate that it becomes of first-order as long as h > 0. 
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FIG. 1. The D/U phase diagram for the anisotropic 2D XY model in the parameter space of 
interchain coupling h and spin-lattice coupling ry. 
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FIG. 2. The order parameter 6* as a function of r] for h = 0,0.05,0.1 and 0.2. 

It is also interesting to notice that the excitation spectrum may be gapless even for finite 
dimerization in presence of the interchain coupling. From ^k, it is easy to find that a gap 
is not fully opened at the Fermi surface until the dimerization S reaches the same value as 
the interchain coupling h. Qualitatively this gives the physical reason for the conditional 
SP transition when /i > as shown by Fig. |1|: The incomplete opening of the gap by the 
dimerization is unfavorable to an adequate reduction for the electronic energy so that the 
gain from it is not always enough to overcome the increase of the elastical energy. Here 
it deserves to mention that the same property of the excitation spectrum is also present 
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for another dimerized state with the pattern corresponding to a wavevector (vr, 0), i.e., 



Uij = (— see Ref. Thus it is reasonable to predict that the similar conditional SP 
transition will occur if this different dimerization pattern is considered. 

We emphasize that the above results were obtained under the 2D JW transformation, 
where the interchain and intrachain exchanges were treated on the same level although the 
phase factor was approximately averaged. One may expect that the effect of interchain 
coupling, i.e., favoring the existence of the U state, could be adequately considered in this 
treatment. Actually, this point has been reflected in Fig. |l], prominently by the relationship 
between rjc and h in the region h —>■ 0: r]c ^ — 1/ln/i. Even an infinitesimal interchain 
coupling may induce a rapid increase of r]c away from zero, i.e., the U state is largely 
stabilized once the interchain coupling is switched on. 

As mentioned before, a natural extension is to consider the case of the Heisenberg model. 
In this case the U state should imply 2;-axis (and rotationally invariant) antiferromagnetic 
long-range order (LRO) [Q. But we believe that the above qualitative results for the XY 
model, especially the rapid increasing property for t]c away from zero, is not changed by the 
extra Ising terms except that the U phase in Fig. |I] is replaced by more precisely called 
Neel phase (i.e., uniform state with antiferromagnetic LRO) More work is necessary 

for justification. 

Before closing we note that recently there are intensive studies on the material CuGeOs 



which show surprising phenomena under doping [^, and theoretically it has been under- 
stood that the interchain coupling plays a key role for these phenomena p5| , p6| . For their 
better understanding our phase diagram given here (for the pure system), which address 
the effect of the interchain coupling, will be helpful. Furthermore it is of great interest to 
extend this work to study the phase diagram for doped quasi-lD SP system. In addition, we 
also notice that the nonadiabatic effect as refered before is important for really modelling 
CuGeOs P,PJ8|J27[| , so inclusion of dynamic phonons beyond this work is another interesting 



topic for future investigation. 

In conclusion, based on the 2D JW transformation we have studied the ground state 
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D/U phase diagram for an anisotropic 2D XY model in adiabatic limit. It is found that the 
spin-lattice coupling rj must exceed some critical value rjc in order to reach the D phase for 
any finite h. The value r]c has the dependence — 1/ln/i on h in the region /i — > and the 
transition between U and D phase is of first-order for at least h > 10~^. 
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for helpful discussions and the support of Visitor Program of MPI-PKS, Dresden. This work 
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